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Abstract 

The warming trend of the last decades is now so strong that it is discernible in 
local temperature observations. This opens the possibility to compare the trend to 
the warming predicted by comprehensive climate models (GCMs), which up to now 
could not be verified directly to observations on a local scale, because the signal- 
to-noise ratio was too low. The observed temperature trend in western Europe over 
the last decades appears much stronger than simulated by state-of-the-art GCMs. 
The difference is very unlikely due to random fluctuations, either in fast weather 
processes or in decadal climate fluctuations. In winter and spring, changes in atmo- 
spheric circulation are important; in spring and summer changes in soil moisture 
and cloud cover. A misrepresentation of the North Atlantic Current affects trends 
along the coast. Many of these processes continue to affect trends in projections for 
the 21st century. This implies that climate predictions for western Europe probably 
underestimate the effects of anthropogenic climate change. 



1 Introduction 

Global warming has been detected in the global mean temperature and on conti- 
nental - scale regions, and th is warming has been attributed to anthropogenic causes 
(|Stottl . I2OO3I: llPCcl l2007l). The ob served global warming trend agrees well with 
predictions (jRahmstorf et al.l . 120071 ). However, climate change projections are typ- 
ically made for much smaller areas. The Netherlands, for instance, corresponds to 
a single grid box in most current climate models, but the temperature projections 
in the KNMFOB scenarios ( van den Hurk et al.L I2OO6II2OO7I) are based on grid point 
values of global and regional climate models. In this region, temperatures simulated 
by Regional Climate models (RCMs) do not deviate much from GCMs, as the pre- 
scribed SST and boundary condition determine the temperature to a large extent 



(jLenderink et al.l . 120071) 



By now, global warming can be detected even on the grid point scale. In 
this paper we investigate the high temperature trends observed in western Europe 
over the last decades. First we compare these with the trends expected on the 
basis of climate model experiments. These turn out to be incompatible with the 
observations over large regions of Europe. The discrepa.ncy is very unlikely due 



to weather or decadal climate fluctuations ( Smith et al.l . 120071 : iKeenlvside et al 



20081 ). Searching for the causes of the unexpectedly fast temperature rise in Europe, 
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we discuss the differences between modelled and observed atmospheric circulation, 
ocean circulation, soil moisture and radiation, aerosols, and snow cover. 



Data 



Many of the results below are obtained in the ESSENCE project, a large ensem- 
ble of climate experiments aimed to objtain a good estimate of internal climate 
variability and extremes |sterl et all l2008h . The ESSENCE database contains re- 
sults of a 17 - member ensemb le of climate runs using the ECHAM5/MPI-0M cli- 
mate model ( Jungclaus et al.l .[2006) of the Max-Planck-Institute for Meteorology in 
Hamburg. The version used here is the same us ed for chmat e scenario runs in prepa- 
ration of the IPC C Fourth Assessment Report (|lPCCl . l2007t ). The ECHAM5 version 
( Roeckner et al 112003) ^las a horizontal resolution of T63 an d 31 vertical hy brid lev- 
els with the top level at 10 hPa. The ocean model MPI-OM (jMarsland et al.l . l2003l ) 
is a primitive equation z-coordinate model. It employs a bipolar orthogonal spher- 
ical coordinate system in which the two poles are moved to Greenland and West 
Antarctica, respectively, to avoid the singularity at the North Pole. The resolution 
is highest, 0(20 — 40 km), in the deep water formation regions of the Labrador, 
Greenland, and Weddell Seas, and along the equator the meridional resolution is 
about 0.5°. There are 40 vertical layers with thickness ranging from 10 m at the 
surface to 600 m at the bottom. 

The experimental period is 1950-2100. For the historical part of this period 
(1950-2000) the concentrations of greenhouse gases (GHG) and sulphate aerosols 
are specified from o bservations, while f or the future part (2001-2100) they follow 
SRES scenario Alb (|Nakicenovic et aLLlioOoh . The runs are initialised from a long 
run in which historical GHG concentrations have been used until 1950. Different 
ensemble members are generated by disturbing the initial state of the atmosphere. 
Gaussian noise with an amplitude of 0.1 K is added to the initial temperature field. 
The initial ocean state is not perturbed. 

The findings from the ESSENCE ensemble are backed with results from en- 
sembles from the World Climate Research Programme's (WCRP) Coupled Model 
Intercomparison Project phase 3 (CMIP3) multi-m odel dataset. We used the mod- 
els wi th the most realistic circulation selected in Ivan Ulden and van Oldenborgh 
( 20061 ). The criterion used was that the explained variance of monthly sea-level 
pressure fields should be positive for all months. The explained variance is given by 



Miff 



(1) 



obs 



Here, cr^jg is the spatial variance of the difference between simulated and observed 
long-term mean pressure, and 0"^,^^ the spatial variance of the observed field. A neg- 
ative explained variance indicates that the monthly mean sea-level pressure deviates 
more from the observed field than the reanalysed field deviates from zero. 

Apart fro m ECHAM5/MPI-0M , the models that were selected are the GFDL 



CM2.1 model ( Delworth et al 



2004 ). HadGEMl (jJohns et al 



2006), MIROC 3.2 T106 (lK-1 model developers 



20041) and CCCMA CGCM 3.2 T63 ( Kim et al 



20021 ). Lower-resolution versions of these models also satisfy the criterion, but were 
thought not to contribute additional information. Observed greenhouse gas and 
aerosol concentrations were used up to 2000, afterwards the SRES Alb scenario 
was prescribed. 
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The model resu lts are compare d with analysed obse rvations in the CRUTEM3 
([Brohan et al.l . l2006l ) and HadSST2 (jRavner et al.l . l2006l) datasets. These have been 
merged with weighing factors proportional to the fraction of land and sea in the grid 
box. For the global mean temperature the variance- weighed HadCRUTS dataset has 
been used. However, this weighing procedure was found to give unrealistic trends 
in the gridded HadCRUTS dataset over Europe in summer. The variance of the 
HadSST2 grid boxes that are mainly land is very small, so these dominate the 
combined value, severely down-weighing the CRUTEM3 land observations. 



3 Trend definition 

Trends are computed as the linear regression against the globally averaged tem- 
perature anomalies, smoothed with a 3yr running mean to remove the effects of 
ENSO, over 1950-200 7. This definition is physically better justified than a linear 
trend (as used in, e.g., Scherrer et al. . 20051 ). and gives a better signal-to- noise ratio. 



In other words, we assume that the local temperature is proportional to the global 
temperature trend plus random weather noise: 

T'{x,y,t) = A{x,y)T'J^ljt) + e{x,y,t) . (2) 

The difference between observed and modelled trends is described by z-values. 
These are derived from the regression estimates and their errors: 

(3) 



y^(AAobs)2 + (A^„,od)V^^ 

with TV the number of ensemble members and the bar denoting the ensemble av- 
erage. The standard errors AA are computed assuming a normal distribution of 
the trends A. The normal approximation has been verified in the model, where the 
skewness of the 17 trend estimates is less than 0.2 in almost all areas where z > 2 
in Fig. [51 Serial correlations have been taken into account whenever significant. 



4 Observed and modelled trends 

Fig. [1^ shows the global annual mean temperature anomalies from observations 
(HadCRUTS) and in the 17-member ESSENCE project ensemble. The model is 
seen to give a very good description of the warming trend so far; the regression of 
modelled on observed global mean temperature is 1.06 ± 0.06. 

In Fig. [TJd the temperature at the model grid point representing the Nether- 
lands is compared with observations at Do B ilt, corrected for change s in observation 



practices and warming due to urbanisation IjBrandsma et al.l . l2003f ). Random fluc- 
tuations due to the weather are much larger at this small spatial scale. In contrast 
to the global trends, the local observations show a much stronger warming trend 
than simulated by this climate model over the last two decades. The model simu- 
lates a factor 1.24±0.09 faster warming than the global mean, but the observations 
have a trend A = 2.23 ± 0.36. The histogram of Fig. [T]; shows that not a single 
ensemble member has a trend this high over 1950-2007. The four other selected 
CMIP3 models also show a trend that is much lower than observed. 
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Figure 1: Annual mean temperature anomalies [K] relative to 1951-1980 in obser- 
vations (red) and the ESSENCE ensemble (blue, 17 realisations and the ensemble 
mean), a) Global mean, b) De Bilt, the Netherlands (52°N, 5°E). c) modelled and 
observed trends [K/K] at De Bilt in the ESSENCE ensemble (histogram) and the 
four other selected CMIP3 climate models. 
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Figure 2: Observed trends in surface temperature (colour, [K/K] March 1950 - 
February 2008, in the merged HadSST2/CRUTEM3 dataset. a) Dec-Feb, b) Mar- 
May, c) Jun-Aug, d) Sep-Nov. A value of one denotes a trend equal to global mean 
warming. The contours indicate the 2; = 2,3 and 4 lines of the significance of the 
difference with the modelled trends (ESSENCE ensemble). Black (red) indicates 
that the observed trend is significantly larger (smaller) than the modelled trend. 



4 



Maps of the observed warming trends A{x, y) in Europe over 1950 to 2007 are 
shown in Fig. [2] As the mechanisms vary over the seasons these are shown separately. 
We also show z-values for the differences between modelled and observed trends by 
contours starting at z = 2. The areas for which Izj > 2 correspond to regions where 
the hypothesis that the model describes the observed trends well can be rejected 
at the 95% confidence level. This area almost coincides with the region where the 
observed trends are higher or lower than any in the 17-member ensemble. 

In all seasons the eastern Atlantic Ocean has warmed significantly faster than 
the model simulated. In spring there are also discrepancies of up to 3 standard 
deviations over land from France to the Baltic and Russia. In summer, the largest 
discrepancies are in the Mediterranean area, the z = 2 contour extending north 
to the Netherlands. In autumn, over land only Great Britain has 95% significant 
discrepancies between observed and modelled trends. 

The area inside the z = 2 contour, 12% to 29% of the area enclosed in 32°- 
72°N, 25°W-35°E, is much larger than the 6% expected by chance at 95% confi- 
dence. For the z = 3 contour the area is 2% to 6%, larger than the 2.5% expected 
except in winter. The area expected by c hance includes the effects of spatial corre- 
lations, assuming 30 degrees of freedom (jLivezev and Chenl . [19831. 

We performed similar analy ses for four other models used for the IPCC Fourth 
Assessment Report (llP Cd, l20^ ) that simulate the current climate in Europe well 



van Ulden and van O ldcnborgh. 20061 ). In Fig. [3] the local temperature trends over 



1950-2007 are shown over Europe in the observations, the ESSENCE ensemble 
of ECHAM5/MPT0M model runs, GFDL CM2.1, MIROC 3.2 T106, HadGEMl 
and CCCMA CGCM 3.2 T63 models. For the models, we define the trend as the 
regression against the modelled global mean temperature0. Over western Europe, 
the patterns of change are similar to the ones in Fig. [21 although the statistical 
significance is lower due to the smaller ensemble sizes. Other ensemble simulations 
appear to exhibit similar behaviour (M. Collins, private communication). Time 
slice experiments of the PRUDENCE ensemble of high-resolution regional climate 
models show temperature cha i iges t hat are similar to the equivalent GCM changes 
( Christensen and Christensenl . l2007h . 

Figs I2|3I show that the probability is very low that the discrepancy between 
observed and modelled warming trends is entirely due to natural variability: the 
area enclosed by the contours is much larger than expected by chance. We therefore 
investigate which physical trends are misrepresented in the GCMs. 



5 Atmospheric circulation 

In Europe, at the edge of a continent, chang es in temperature are caused to a large 
extent by changes in atmospheric circu l atioii ( Osborn and Jone^ 2000t Turnpennv et al 



20021 : Ivan Oldenborgh and van uidi3,[20o3). To investigate the effects of trends in 



the atmospheric circulation, monthly mean temperature anomalies are approxi- 
mated by a simple model that i solates the linear effec t of ci rculation anomalies 
( van Ulden and van Oldenborg E 120061 : I van Ulden et al 1 l2007h . These are the ef- 



fects of the mean geostrophic wind anomalies U'{t),V'{t) across the temperature 
gradients, and vorticity anomalies W'{t) that influence cloud cover. The other terms 
are the direct effect of global warming, approximated again by a linear dependence 



1. The MIROC 3.2 T106, HadGEMl and CCCMA CGCM 3.2 T63 experiments in the CMIP3 
archive exhibit an 0(1.5) times faster global mean temperature rise than observed. 
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Figure 3: The trends in temperature in western Europe as the regression against 
(modened) global mean temperature [K/K] in the observations and the GCMs with 
the most realistic mean circulation in Europe over 1950-2007. The contours denote 
the number of standard errors between the observed and modelled trends starting 
at z = 2 (black) and z = —2 (red). 
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Figure 4: As Fig. [21 but for the circulation-dependent temperature Tci 



on the global mean temperature TgiobailOi the remaining noise rjit). A memory 
term M describes the dependence on the temperature one m onth earlier, which is 
important near coasts ( van Ulden and van Oldenborgh . [20061 ): 
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The geostrophic wind and vorticity anom alies U',V\W' are computed from the 
NCEP/NCAR reanalysis sea-level pressure ( Kalnav et al. . 1996[ ) and the coefficients 
M, Ay, Ajj, B and A are fitted over 1948-2007 for each calendar month. This model 
explains more than half the variance in monthly mean temperature over most of 
Europe, both in the observations and the models (with coefficients fitted from model 
data). Temperature changes that are due to changes in the atmospheric circulation 
show up as trends in T^'j^^.. Fig. 3] shows the warming trends in the circulation- 
dependent temperature in the observations and the significance of the difference 
with the ECHAM5/MPI-0M chmate model results. 

In winter, the observed temperature rise around 52°N is dominated by cir- 
culation change s. Fig. [5^ show s that a significant increase in air pressure over the 
Mediterranean ( Osborrl 20041 ) {z > 3) and a not statistically significant air pres- 
sure decrease over Scandinavia (z < 2) have brought more mild maritime air into 
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Europe north of the Alps. 

In Fig. [5] trends in sea-level pressure over 1950-2007 of the NCEP/NCAR 
reanalysis are compared to climate model simulations. Both the reanalysis and the 
ESSENCE ensemble show a significant trend in the Mediterranean region, but the 
observed trend is a factor four larger than the modelled trend. The GFDL CM2.1 
and MIROC 3.2 T106 models also show significant positive trends in this area, but 
again much smaller than observed. The other two models show no positive trends 
there. 

We conclude that the temperature trends in winter and to a lesser extend 
spring are due to a shift towards a more westerly circulation. This change is under- 
represented in climate models. In summer and autumn the rise in temperature is 
mainly caused by factors not linearly related to shifts in atmospheric circulation. 



6 Oceanic circulation 

The temperature trend in the eastern Atlantic Ocean is underestimated by the 
model results in all seasons but summer and this motivated an investigation of the 
Atlantic ocean circulation. The discrepancy may be either a result of ocean memory 
of the initial state, or model errors. 

The ESSENCE ensemble was started from a common ocean initial state in 
the model year 1950. This initial state was taken from a coupled run, so it does 
not correspond to the real state of the ocean in 1950. It has recently been shown 
that ocean mem ory and dynamics lead to predictability in year 5-10 in the North 
Atlantic Ocean ( Keenlvside et al. . 20081) . However, after 10 years the ocean states 



have decorrelated completely, as is illustrated by the autocorrelation function of 
the maximum overturning circulation at 35°N and an index of the Atlantic Multi- 
decadal Oscillation (AMO) shown in Fig. \6\ This result is in agreement with the 
decorr elation time of less than 10 year s found in a large ensemble of the CCSM 1.4 
model ( Driifhout and Hazeleger . 20071 ). As our definition of the trends does not give 



weight to temperature variations in the first ten years, when the global mean tem- 
perature is almost constant, the effect of ocean memory on the trends is negligible. 
The fact that the observed trend is outside the ensemble spread therefore includes 
the effects of decadal climate variations, to the extent that these are simulated well 
by the models. 

In the observations the multi-decadal oscillations in the Atlantic Ocean are 
stronger and slower (Fig. [6]) than in the ECHAM5/MPI-0M model. Over the last 
decades there has been a rising trend in the AMO index. To disentangle the effects 
of the AMO and global warming on temperatures in the North Atlantic region, we 
subtract a term proportional to the global mean temperature from the SST average, 
fitted over the 150 years with estimates for both. In the model, this gives the same 
result as subtracting the ensemble mean (the AMO has very little effect on the 
global mean temperature). Over the relatively short period 1950-2007 we then find 
virtually no contribution from the AMO on the trend in the observations either. 

Systematic model errors play a much larger role. The coarse resolution ocean 
models used in GCMs have a common error in the North Atlantic Current (NAC). 
The NAC is compa red between the 0.5° SODA 1.4.1 and 1.4.2 ocean reanalyses 
(ICarton et al.l . l2005l) and the ECHAM5/MPI-0M GCM. Fig.[7t,d show that in the 



average over the upper 750m, the warm water of the modelled NAC crosses the basin 
zonally to Portugal, and continues northward, whereas in the reanalysis this Azores 
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Figure 5: Trends in Dec-Feb sea-level pressure [hPa/K] over 1950-2007 in the 
NCEP/NCAR reanalysis (a), ECHAM5/MPI-0M (b), GFDL CM 2.1 (c), MIROC 
3.2 T106 (d), CCCMA CGCM 3.1 T63 (e) and HadGEMl (f). The contours denote 
the z-value of the trend being different from zero, starting at 2. 
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Figure 6: Autocorrelation function of the ECHAM5/MPI-0M Atlantic Meridional 
Overturning Circulation (AMOC) at 35°N and the Atlantic Multidecadal Oscilla- 
tion (AMO) index, SST averaged over 25°-60°N, 75°-7°W. The effects of external 
forcing have been minimised by taking anomalies relative to the ensemble mean in 
the model, and by subtracting the regression against the global mean temperature 
in the observations. 



current is much weaker and most water meanders north-east across the Atlantic 
as part of the surface bran ch of the Atlantic Meridional Overturning Circulation 
(jLumpkin and Speeii . 1200,'^ . 

The mean vertical thermal structure is shown in Fig. [7l3,f at 55°N. The bias 
in the currents results in a too weak vertical stratification and very deep mixed 
layers in the modelled East Atlantic, where the surface is cooled by cold fresh 
water advected from the north (due to too strong westerlies that drive a too large 
southward Ekman drift, Figs[7K,b) and warmed by the anomalously warm water 
below (associated with a too far eastward flowing NAC). The deep mixed layer 
hardly warms under global warming, whereas observed surface temperature rises at 
about the same rate as the global mean temperature. 

A signature of this bias in the NAC is a strong negative SST bias in the 
middle of the northern Atlantic Ocean. In the observations this region is south of 
the NAC, but in the models it is located north of the current and hence it is much 
colder. Such a bias is clearly visible in all CMIP3 models considered (Fig.[8|3-f), but 
absent when comparing a the high-res olution SODA reanaly sis to the same lower 
resolution Oi v2 SST analysis fFig. [8k iRevnolds et all . l2002f ). 

This bias in the ocean explains the discrepancies over the ocean in Figs I2|3I 
To estimate the effect on land temperatures, we approximated the effect of a bias in 
the trend in the East Atlantic on 2-m temperature in Europe by the effect of decadal 
variability in the same region in the ESSENCE ensemble over 1950-2000. For each 
month, the regression of 2-m temperature was computed with SST averaged over 
40°-50°N, 30°-10°W the previous month, low-pass filtered with a Syr running mean. 
Trends were removed by taking anomalies with respect to the 17- member ensemble 
mean. The results are shown in Fig. [9l There is an influence of East Atlantic SST 
on coastal temperatures of 0.3 to 0.5 K per degree change of East Atlantic SST the 
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Figure 7: Ocean surface currents [ms^^] in the SODA reanalysis (a) and the 
ESSENCE ensemble mean (b), both averaged over 1961-1990. Northward currents 
are shown positive, southward currents negative, the colour denotes the total ve- 
locity. The same for vertically integrated currents from to 750m [m^s^^] (c,d). 
Subsurface temperature [°C] across the Atlantic Ocean at 55°N in SODA (e) and 
the ESSENCE ensemble (f). 
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Figure 8: Difference between 1982-2007 annual mean SST and the 01 v2 SST analy- 
sis: SODA ocean reanalysis (a), ESSENCE ensemble (b), GFDL CM 2.1 (c), MIROC 
3.2 T106 (d), HadGEMl (e) and CCCMA CGCM 3.1 T63 (f). 




Figure 9: Regression of local temperature on SST averaged over 40°-50°N, 30°- 
10°W in the ESSENCE ensemble, low-pass filtered with a 5yr running mean, sum 
of monthly 1-month lag regressions with SST leading, 1950-2000, anomalies w.r.t. 
the ensemble mean. Dec-Feb (a), Jun-Aug (b). 
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previous month, but the signal does not extend very far inland. 



7 Soil moisture and short-wave radiation 

The third important factor explaining discrepancies between observed and modelled 
trends in Figs I2I3I consists of related trends in soil moisture and radiation at the 
surface in spring and summer. In summer, the pattern of stronger-than-expected 
heating corresponds closely to the area in which evapotrans p iratio n correlates neg- 
atively with temperature in the RCM of 'S eneviratne et al. I (I2OO6I ) (their Fig. 3a). 



This indicates that in this area, the soil moisture is exhausted to the extent that 
an increase in radiation translates directly into a large increase in temperature, 
whereas in wetter areas the evapotranspiration increases with rising temperature, 
damping the high temperatures. It should be noted, however, that the observed 
trend (2.6 ±0.2 over 40°-50°N, 0°-15°E) is much stronger than the modelled trend 
(1.4 ±0.1), indicating that the GCMs underestimate the strength of this process in 
the current climate. 

Regional climate models do not resolve t his discrepancy. Comparing the ES- 



SENCE results with the PRUDENCE ensemble (IChristensen and Christensenl . 120071 ). 
we find that the second-highest temperature increases in the Mediterranean, the 
Alps and southern France between 1960-1990 and 2071-2100 are no more than 
25% higher than the equivalent numbers for ECHAM5/MPI-0M, whereas the dis- 
crepancy between observed and modelled trends approaches a factor two. There 
is therefore no indication that RCMs simulating the last 50 years would show a 
warming trend as high as observed. 

To explain the warming trends further n orth, we propos e a me chanism that 
closely resembles the mechanism described in Vautard et al.l ( 2007 ) for extreme 



summers in Europe. North of the area with most severe drying, southerly winds 
bring warmer and drier air northwards, increasing the amount of solar radiation 
reaching the ground. Northerly winds do not change. With the wind direction ran- 
domly fluctuating between these two, the net effect is a heating trend accompanied 
by soil drying. This way the effects of soil moisture depletion migrate northwards. 

We found supporting evidence using Dutch global short-wave radiation obser - 
vations, which are well-calibrated since the early 1970s (iFrantzen and Raafi . 1978h . 



The monthly mean observations were corrected for circulation effects using a model 
analogous to Eqs Q-®. The trend in circulation is small in late spring and summer 
(cf. Fig.|4]), so subtracting circulation effects mainly decreases the variability. 

All six stations with observations show an increase in global short-wave radia- 
tion in spring and summer over the period 1971-2007, averaging to 14±2Wm^^K^^ 
(Fis. [T0|). To translate changes in short-wave radiation to temperature changes we 
use a conversion factor obtained from the regression of detrended monthly mean 
temperature on incoming short-wave radiation, which is 0.05 K/Wm~^. The ob- 
served long-term trend in global short-wave radiation corresponds to roughly 0.7 K 
warming per degree global mean temperature rise. This is a sizeable fraction of the 
total temperature trend, 3.0 ± 0.5 K/K in spring and 2.2 ± 0.6 K/K in summer. 

The GCM also has a positive trend in this area, but only 5±2Wm-2K-i over 
1971-2007. The difference, equivalent to a trend of 0.5 in units of global mean tem- 
perature, therefore explains half the discrepancy between observations and model 
in the Netherlands. Spatially, the modelled trend in short-wave radiation is at the 
northern side of the area of strongest warming in Fig. [2t, in accordance with our 
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Figure 10: Trends over 1971-2007 in global short-wave radiation [Wm ^] in 
spring (a) and summer (b) in the ESSENCE ensemble of 17 ECHAM5/MPI-0M 
model experiments, the Netherlands average, and all stations in the Netherlands. 
Error bars denote the standard error. 



hypothesis for the summer. In the model the trend is mainly due to a decrease 
in cloud cover and continues up to 2100, also supporting the hypothesis that the 
decrease in cloudiness is driven by soil moisture depletion further south. We do not 
have an explanation for the increased sunshine in spring. 

There are indications in the observations that the trend is largest on days 
with southerly wind directions, both in spring and in summer, but the statistical 
uncertainty on these res ults is large. Direct cloud cover observations are unreliable 
(iNorris and Wildl. l2007^ and uncertainties in cloud cover changes are known to be 
large (llPCd . l2007ir making this mechanism difficult to investigate further using 
observations, but likely to be relevant. 

Land use changes are estimated to contribute 0(0. IK) to the temperature rise 
in the Netherlands up to now. This value comes from a dir ect estimate of the effect 
of growing cities around De Bilt ( Brandsma et al.l . l2003f l. A rough country- wide 
estimate can be deduced from t he measured increase in 'built-up area' of 1%/lOyr 
over 1986-1996 and 1996-2003 (jCentraal Bureau voor Statistied . l2007l ). Assuming 
the latent heat flux is halved over this area, this decreases evaporative cooling 
by 0{2 Wm~^) over 30 years, causing a 0(0.1 K) temperature rise. We conclude 
that land use changes do not contribute substantially to the discrepancy between 
observed and modelled temperature trends. 



8 Aerosols 



Air pollution has decreased summer tem peratures in Europe from 1950 to around 
1985 , after this clearer skies (ISternl. 2006h have c aused a temperature rise ( Wild et al 



20051 : iNorris and Wildl . l2007t IWild et al.l . l2007f ). This is reflected in first a decrease 



and later an increase in observed short-wave radiation of about 0.3 Wm^ yr~ in 
the Netherlands in summer (see Fig. Ilip . Converting to an annual mean, this is on 
the low end of the range quoted European average of 0.3±0.1 Wm~^yr~^, corrected 
for cloud cover changes (Norris and Wild, 2007). As the Netherlands, on the coast, 
escaped the worst affects of air pollution, this difference is not surprising. 
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Figure 11: Modelled global circulation- independent short-wave radiation [Wm~^] 
compared with observations at the two stations with the longest records in the 
Netherlands in spring (a) and summer (b). 



The observed decrease over 1970-1985 translates into a cooling effect of 0.3 to 
0.4 K. Note that the effect of this temporary dimming on the trend over the longer 
period 1971-2007 or 1950-2007 is small: the dimming and brightening cancel each 
other to a large extend. 

In our trend measure the effect of decreased solar radiation due to direct 
and indirect aerosol effects is about 0.2 times the global mean temperature. This 
explains only a small part of the observed trend in the Netherlands in summer. On 
shorter time scales, e.g. the period 1985-2007, the reduction of aerosols of course 
gives a much larger contribution to the temperature trend. 

The incoming solar radiation in the ESSENCE ensemble shows a smaller 
aerosol effect of 0.1 ± 0.1 Wm~^yr~^ in the Netherlands in summer. The discrep- 
ancy translates into a temperature trend bias of only 0.1 ± 0.1 K per degree global 
warming, significantly smaller than the effect of the bias in long-term trend dis- 
cussed above. 
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9 Snow cover 

In spring, differences in modelled and observed snow cover trends amplify the dis- 
crepancies in trends in the Baltic region. In Fig. [12] the trend in Mar-May snow 
cover is shown in the observations and the ESSENCE ensemble. The observations 
indicate a much faster decrease of spring snow cover than the model. At most grid 
points the significance of the difference is not very high (p < 0.2) because of the 
large decadal fluctuations in the observed snow cover. 



10 Conclusions 

We have shown that the discrepancy between the observed temperature rise in 
western Europe and the trend simulated in present climate models is very unlikely 
due to fast weather fluctuations or decadal climate fluctuations. The main physical 
mechanisms are varied, both geographically and as a function of the seasonal cycle. 
The most important discrepancies between observations and models are 

1. a stronger trend to westerly circulation in later winter and early spring in the 
observations than in the models, 

2. a misrepresentation of the North Atlantic Current in the models giving rise 
to an underestimation of the trend in coastal areas all year, and 

3. in summer, higher observed than modelled trends in areas in southern Europe 
where soil moisture depletion is important, and consequently a stronger trend 
in sunshine around the Netherlands in spring and summer. 

Smaller contributions come from differences between observed and modelled trends 
in aerosol effects in spring and summer, and snow cover changes in the Baltic in 
spring. 

As most projections of temperature changes in Europe over the next century 
are based on GCMs and RCMs with the biases discussed above, these projections 
are probably biased low. To correct the biases, it is essential to not only validate 
the GCMs for a good representation of the mean climate, but also on the observed 
temperature trends at regional scales. 
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